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Abstract 

In this paper we provide sharp criteria for hnear stabihty or instabihty of equihbria of 
colhsionless plasmas in the presence of boundaries. Specifically, we consider the relativistic 
Vlasov-Maxwell system with specular reflection at the boundary for the particles and with the 
perfectly conducting boundary condition for the electromagnetic field. Here we initiate our 
investigation in the simple geometry of radial and longitudinal symmetry. 

1 Introduction 

We consider a plasma at high temperature or of low density such that collisions can be ignored as 
compared with the electromagnetic forces. Such a plasma is modeled by the relativistic Vlasov- 
Maxwell system (RVM) 

f dtf+ + V ■ + {E + vxB)- = 0, 

[ dtr + V ■ VJ- - (E + * X B) • V,/- = 0, 

V^-B = p, V^-B = 0, (1.2) 
9tE - X B = -j, a^B + X E = 0, (1.3) 

{f+-r)dv, j= / v{f+-r)dv. 



Here f^{t,x,v) > is the density distribution for ions and electrons, respectively, x € O C 
is the particle position, O is the region occupied by the plasma, v is the particle momentum, 
(v) = \/l + |f p is the particle energy, v = v/{v) the particle velocity, p the charge density, j the 
current density, E the electric field, B the magnetic field and ib(E + v x B) the electromagnetic 
force. We assume that the particle molecules interact with each other only through their own 
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electromagnetic forces. For simplicity, we have taken all physical constants such as the speed of 
light and the mass of the electrons and ions equal to 1. This whole paper can be easily modified to 
apply with the true physical constants. 

Stability analysis for a Vlasov-Maxwell system of the type that we present in this paper has 
so far appeared only in the absence of spatial boundaries, that is, either in all space or in a 
periodic setting like the torus. In this paper we present the first systematic stability analysis in a 
domain il. with a boundary. It is an unresolved problem to determine which boundary conditions 
an actual plasma may satisfy under various physical conditions. Several boundary conditions are 
mathematically valid and some of them are more physically justified than others. Stability analysis 
is a central issue in the theory of plasmas. In a tokamak and other nuclear fusion reactors, for 
instance, the plasma is confined by a strong magnetic field. This paper is a first, rather primitive, 
step in the direction of mathematically understanding a confined plasma. We take the case of a 
fixed boundary with specular and perfect conductor boundary conditions in a longitudinal and 
radial setting. 

The specular condition is 

f^{t, X, v) = f^{t, x,v - 2{v ■ n{x))n{x)), n{x) • i; < 0, x G 50, (1.4) 

where n{x) denotes the outward normal vector of 00 at x. The perfect conductor boundary 
condition is 

Exn(a;) = 0, B • = 0, x S (90. (1.5) 
Under these conditions it is straightforward to see that the total energy 

£it) = lf [ v\f+ + r)dvdx + l [ (\B\'' + \B\^) dx (1.6) 

is conserved in time, and also that the system admits infinitely many equilibria. The main focus 
of the present paper is to investigate stability properties of the equilibria. 

Our analysis closely follows the spectral analysis approach in |15^ [T6\ [T7] which tackled the 
stability problem in domains without spatial boundaries. Roughly speaking, that approach provided 
the sharp stability criterion > 0, where is a certain nonlocal self-adjoint operator acting on 
scalar functions that depend only on the spatial variables. The positivity condition was verified 
explicitly for various interesting examples. It may also be amenable to numerical verification. In 
our case with a boundary, every integration by parts brings in boundary terms. This leads to some 
significant complications. 

In the present paper, we restrict ourself to the stability problem in the simple setting of longi- 
tudinal and radial symmetry. Thus the problem becomes spatially two-dimensional. Indeed, using 
standard cylindrical coordinates {r,9,z), the symmetry means that there is no dependence on z 
and 9 and that the domain is a cylinder O = D x M where D is a disk. We may as well assume that 
D is the unit disk in the (xi, X2) plane. It follows that x = {xi,X2, 0) S B x M, t; = (i;i, ^2, 0) G M^, 
E = {El, E2,0), and B = (0, 0, B). In the sequel we will drop the zero coordinates so that x € B and 
V G M^. In terms of the polar coordinates {r,6), we denote = {cos6,sm9), eg = (— sin 0, cos 0). 
It follows that the field has the form 

E = -dr(per - dtil^eg, B = -driril))), (1.7) 
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where the scalar potentials <p{t, r) and "0(4, r) satisfy a reduced form of the Maxwell equations. See 
the next section for details. 



1.1 Equilibria 

We will denote an equilibrium by (/^'^, E*^, Its field has the form 

r 

Then the particle energy and angular momentum 

e^(x, v) := (v) ± (p^{r), p^{x, v) := r{vg ± 'ip'^{r)), 



(1.9) 



are invariant along the particle trajectories. It is straightforward to check that fi^{e^,p^) solve 
the Vlasov equations for any smooth functions fi^{e,p). So we consider equilibria of the form 

= fi~^{e'^{x,v),p'^{x,v)), f~{x,v) = fi" {e~ {x , v) , p~ {x , v)) . (1.10) 

The potentials still have to satisfy the Maxwell equations, which take the form 



-A/ 

-A^V° 



/i+(e+,p+) - fj. {e ,p ) 



dv 



Vg 



/i+(e+,p+) - {e ,p ) 



(1.11) 



dv 



with Ar = A — ^. Again, see the next section for details. It is clear that the boundary conditions 
(jl.4p and (jl.Sp are automatically satisfied for the equilibria since and p^ are even in Vr, and E*^ 
is parallel to er- In the appendix we will prove that plenty of such equilibria do exist. 

Let (Z'^'^, E", B*^) be an equilibrium as just described with f^'^ = fi^{e^,p^). We assume that 
/i^(e,p) are nonnegative, smooth, and satisfy 



fif{e,p)<0, \fi^{e,p)\ + \fi^{e,p)\ + 



\l^e{e,p)\ 



< 



C„ 



1 + ep 



(1.12) 



for some constant and some 7 > 2, where the subscripts e and p denote the partial derivatives. 
In addition, we also assume that ip^,ip^ are continuous in D. It follows that E'^jB'^ € C^(B), as 
proven in the appendix. 

We consider the Vlasov-Maxwell system linearized around the equilibrium. The linearization is 



dtf^ + D±/± = ^(E + {) X B) • V,/0'±, 



(1.13) 



together with the Maxwell equations and the specular and perfect conductor boundary conditions. 
Here denotes the first-order linear differential operator: := v ■ it (E" + vx B*^) • V^,. See 
the next section for details. 
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1.2 Spaces and operators 

In order to state precise results, we have to define certain spaces and operators. We denote by 
n± = Lf ± (D X ]R2) the 

weighted space consisting of functions f^{x,v) which are radially 

Me I 

symmetric in x such that 

IMc^II/^P dvdx < +00. 



The main purpose of the weight function is to control the growth of as |u| ^ 00. Note that the 
weight ln^l never vanishes and it decays like a power of v as l^l — ?• 00. When there is no danger of 
confusion, we will often write T-L = 7i^. 

For A: > we denote by HI^CB) the usual space on D that consists of functions that are 
radially symmetric. If A; = we write L^(D). By H^^{n) we denote the space of functions = ilj{r) 
in i7^(D) such that il){r)e^^ belongs to the usual H^iJ])) space. The motivation for this space is to 
get rid of the apparent singularity 1/r^ at the origin in the operator A^, thanks to the identity 

By V we denote the space consisting of the functions in i?^(B) which satisfy the Neumann boundary 
condition and which have zero average over D. Also, let V+ be the space consisting of the functions 
in ff^^(D) which satisfy the Dirichlet boundary condition. The spaces V and incorporate the 
boundary conditions (j2.9p for electric and magnetic potentials, respectively. 

Denote by the orthogonal projection on the kernel of in the weighted space n^. In the 
spirit of [151 117j . our main results involve the three linear operators on L^(B), two of which are 
unbounded, 

A\h = -Ah - j ^+(1 - V+)h dv- j r~)h dv, 

A^h = -Arh- j rveifi^ + Hp) dvh - j vg(^fxtV^ (veh) + fJ-^V" (veh)^ dv, (1.14) 
B°h = r / + dvh+ / (fit'P'^iveh) + n~V~{veh)) dv. 



Here denotes fi^{e^,p^) = fJ-^{{v) ± ip^{x),rve ± rilj^{x)). These three operators are naturally 
derived from the Maxwell equations when /+ and /~ are written in integral form by integrating 
the Vlasov equations along the trajectories. See Section [3.21 for their properties. In the next section 
we will show that both Ai with domain V and A2 with domain are self-adjoint operators on 
L^(B). Furthermore, the inverse of Ai is well-defined on the range of B^, and so we are able to 
introduce our key operator 

C° = A'i + {B°)*{A^^r^B^. (1.15) 

The operator will then be self-adjoint on L^(D) with its domain V^. As the next theorem states, 
£0 > [which means that {C^ip, 'iP)l^ ^ for all Tp € V^] is the condition for stability. 

Finally, by a growing mode we mean a solution of the linearized system (including the boundary 
conditions) of the form (e^*/=^, e^*E, e^*B) with ReX > such that G T-L"^ and E,B G L'^{B). 
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The derivatives and the boundary conditions are considered in the weak sense, which will be justified 
in Lemma l2.2i In particular, the weak meaning of the specular condition on will be given by 

1.3 Main results 

The first main result of our paper gives a necessary and sufficient condition for stability in the 
spectral sense. 

Theorem 1.1. Let (Z'^'^, E'', B") be an equilibrium of the Vlasov- Maxwell system satisfying (jl.l2p . 
Consider the linearization of the Vlasov-Maxwell system (|1.13|) for radially symmetric perturbations 
together with the specular and perfect conductor boundary conditions. Then 
(i) if >(), there exists no growing mode of the linearized system; 

(a) any growing mode, if it does exist, must be purely growing; that is, the unstable exponent A 
must be real; 

(Hi) if ^ 0, there exists a growing mode. 

Our second main result provides explicit examples for which the stability condition does or does 
not hold. For more precise statements of this result, see Section [5l 

Theorem 1.2. Let (/i^, E", B") be an equilibrium as above. 

(i) The condition pfJ.^{e,p) < for all {e,p) implies > 0, provided that ip^ is bounded and 
ip^ is sufficiently small. (So such an equilibrium is stable.) 

(a) The condition \fx^{e,p)\ < i^''^^^-, for some 7 > 2 and for e sufficiently small implies > 0, 
provided that ip^ = 0. Here is not necessarily small. (So such an equilibrium is .stable.) 

(Hi) The conditions fi~^{e,p) = ^i~{e,—p) and pf^i~ {e,p) > cop^i^(e), for some nontrivial non- 
negative function iy{e), imply that for a suitably scaled version of (^^,0,B'^), > is violated. 
(So such an equilibrium is unstable.) 

The instabilities in a plasma are due to the collective behavior of all the particles. For a ho- 
mogeneous equilibrium (without an electromagnetic field) Penrose [19] found a beautiful necessary 
and sufficient condition for stability of the Vlasov- Poisson system (VP). For a BGK mode, the equi- 
librium has an electric field. In that case, proofs of instability for VP (electric perturbations) were 
first given in [71 El [9] and especially in |13l [T3] for non-perturbative electric fields. Once magnetic 
effects are included, even for a homogeneous plasma, the situation becomes much more complicated: 
see for instance [6l \TU[ [TT] . In a series of three papers \15\ [TBI HZ] ^ more general approach was 
taken that treated fully electromagnetic equilibria with electromagnetic perturbations. The linear 
stability theory was addressed in |15[ [T7j , while in [16] fully nonlinear stability and instability was 
proven in some cases. In all of the work just mentioned, boundary behavior was not addressed; the 
spatial domains were either all space or periodic. 

In this paper we do not address the question of well-posedness of the initial-value problem. 
For VP in well-posedness was proven in [20] and [T^. For RVM in all space ffi^ it is a famous 
open problem. The relativistic setting seems to be required, for otherwise the Vlasov and Maxwell 
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characteristics would collide. However, for RVM in the whole plane M?, which is the case most 
relevant to this paper, well-posedness and regularity were proven in [3]. The same is true even in 
the 2.5 dimensional case [3]. Although global weak solutions exist in M^, they are not known to be 
unique [2]. Furthermore, in a spatial domain with a boundary on which one assumes specular and 
perfect conductor conditions, global weak solutions also exist [S]. Well-posedness and regularity of 
VP in a convex bounded domain with the specular condition was recently proven in |12j . 

A delicate part of our analysis is how to deal with the specular boundary condition within the 
context of weak solutions. This is discussed in subsection 12. 3i Properly formulated, the operators 

then are skew-adjoint. In this paper, as distinguished from [15], we entirely deal with the weak 
formulation. Our regularity assumptions are essentially optimal. In subsection 12.41 we prove that 
the densities of a growing mode of the linearized system decay at a certain rate as ["yl — t- oo. 

As in [15], the stability part of Theorem 1 1.1 1 is based on realizing the temporal invariants of the 
linearized system. They have to be delicately calculated due to the weak form of the boundary 
conditions. This is done in subsection 13.31 The invariants are the generalized energy X and the 
casimirs ICg, which are a consequence of the assumed symmetry of the system. A key part of the 
stability proof is to minimize the energy with the magnetic potential being held fixed (see subsection 
13. 4p . The purity of any growing mode, in part (ii) of Theorem ll.l|, is a consequence of splitting 
the densities into even and odd parts relative to the variable Vr (see subsection 13. 5p . The proof of 
stability is completed in subsection 13.61 

The proof of instability in Theorem requires the introduction of a family of linear 

operators which formally reduce to as A — ?• 0. The technique was first introduced by 
Lin in [13] for the BGK modes. These operators explicitly use the particle paths (trajectories) 
{X^[s]x,v),V^{s;x,v)). The trajectories reflect specularly a countable number of times at the 
boundary. We use them to represent the densities in integral form, like a Duhamel representa- 
tion. This representation together with the Maxwell equations leads to the family of operators in 
subsection 14.31 Self-adjointness requires careful consideration of the trajectories. It is then shown 
in subsection 14.41 that is a positive operator for large A, while it has a negative eigenvalue for 
A small because of the hypothesis ^ 0. Therefore has a nontrivial kernel for some A > 
If is in the kernel, it is the magnetic potential of the growing mode. From we construct the 
corresponding electric potential ip and densities 

In Section [5] we prove Theorem 11.21 The stability examples are relatively easy. As we see, the 
basic stabilizing condition is pn^ < 0. To construct the unstable examples we make the simplifying 
assumption that the equilibrium has no electric field so that = In the expression [C^ip, '4^)l'^ 
the term that has to dominate negatively is the one with p^p. In order to make it dominate, we 
scale the equilibrium appropriately. We first treat the homogeneous case (Theorem 15. 3p and then 
the purely magnetic case (Theorem 15. 4p . 
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2 The symmetric system 



2.1 The potentials 

It is convenient when dealing with the Maxweh equations to introduce the electric scalar potential 
and magnetic vector potential A through 

E = -Vv?-(9tA, B = VxA, (2.1) 

in which without loss of generality, we impose the Coulomb gauge constraint V • A = 0. Note that 
with these forms, there automatically hold the two Maxwell equations: 

+ V X E = 0, V • B = 0, 

whereas the remaining two Maxwell equations become 

-Av3 = /3, St^A- AA + 5jV(^=j. 

Under the assumption of radial and longitudinal symmetry, there is no z or dependence. We 
use the polar coordinates x = {r cos 6 , r sin 9) on the unit disk B). A radial function f{x) is one 
that depends only on r and in this case we often abuse notation by writing f{r). We also denote 
the unit vectors by er = (cos 0, sin 0) and eg = {— sm6,cos9). Thus er{x) = n{x) is the outward 
normal vector at x € OD. Although the functions do not depend on 0, the unit vectors er and eg 
do. Then we may write = /^(i, r, u), where v = Vrer + vgeg and A = A^e^ + Ageg. 

Now the Coulomb gauge in this symmetric setting reduces to ^^^('^^r-) = 0, so that Ar = h{t) jr. 
We require the field to have finite energy, meaning that E, B G L^(D). Thus Er = —dr^p — dtAr G 
-L^(B) only if h{t) is a constant. But if h{t) is a constant, we may as well choose it to be zero 
because it will not contribute to either E or B. For notational convenience let us write ^ in place 
of Ag. The fields defined through ()2.ip then take the form 

E = -dr^er - dtipeg, B = Be^, B = -9^(r^)). (2.2) 

r 

We note that ^ ^ 

V ■ Vxf = Vrdrf + -Vgdgf = Vrdrf + -Vg^Vgd^J - Vrd^J^. 

In these coordinates the RVM system takes the form 

+ Vrdrf + + (Er + VgB + lvgVg)d,J+ + (Eq - VrB - ^VrVgY^J^ = 0, 
dtf~ + Vrdrf" - [Er + VgB - IvgvgjdyJ'' - (^Eg - VrB + ^VrVg^d^J' = 0, 



(2.3) 



-Aip = p= [ {f+ - f-){t,r,v) dv, 

dtdr^p = jr = / Vrif^ - f~){t, r, v) dv, (2.4) 
(9| - Ar)^ = jg= [ Oe(/+ - f-){t, r, v) dv, 
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where A,. = A — The system (j2.3p - (|2.4p is accompanied by the specular boundary condition 
on which is now equivalent to the evenness of in Uj. at r = 1. In particular, jV = on dD. 
The condition B • n = is automatic. The boundary conditions on if and tp are 

drfit, 1) = const., ip{t, 1) = const. (2-5) 

The Neumann condition on (p comes naturally from the second Maxwell equation in ()2.4p with 
jr = 0. The Dirichlet condition on ip comes naturally from = E x n = (0, 0, —dtip). 

2.2 Linearization 

We linearize the Vlasov-Maxwell near the equilibrium E°, 5°). Fr om (12.31) . the linearized 

Vlasov equations are 

+ D+/+ = -{Er + veB)d,J^'+ - {Eg - VrB)d,,f 
dtr + D-/- = {Er + veB)d,j''-~ + {Eg - VrB)d.,J' 
The first-order differential operators := v ■ ± (E^ + i) x B*^) • V^; now take the form 



0.- , .0,- (2-6) 



D+ := vrdr + + veB° + ^veve)dv^ - (vrB^ + \rVe)d^g, 

D~ := Vrdr - {Er + ^dB^ - ^VgVg^dy^ + (vrB^ - ^VrVg'jdvg. 



(2.7) 



In order to compute the right-hand sides of ()2.6p more explicitly, we differentiate the definition 

Thus, together with the forms of E and B in (j2.2p . we calculate 

{Er + veB)d,,.f'+ + {Eg - VrB)d,,f'+ 

= {Er + VgB)^l^Vr + {Eg - VrB){^4ve + r^i^) 
= -i4vrdrf - HpVrdr{rtp) - fJ-tvedtip - rii^dti) 
= -M^D+(^ - ^+D+(r^) - (^i+ue + rfi+)dtij, 

where the last line is due to the fact that = Vrdrf for radial functions. Of course a similar 

calculation holds for /'^'~. Thus the linearization (12. 6p becomes 

+ D+/+ = /i+D+v? + /x+D+(rV) + 9,,, [^+19^^' 

= fl'^Vrdr^p + HpVrdr{rip) + {fj-tvg + rfj.^)dtip ^^^^ 
9t/" + D~/" = -fi^B^if - fipB~{r'ip) -dyg[n~]dtiJ 

= -fi~Vrdr(p - fl~Vrdr{r'llj) - {^i~vg + rfip)dttp 



Of course, linearization does not alter the Maxwell equations (j2.4p . As for boundary conditions, 
we naturally take the specular condition on and 

dr'p{t,l)=0, 1) = 0. (2.9) 
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2.3 The Vlasov operators 

The Vlasov operators are formally given by (|2.7|) . Their relationship to the boundary condition 
is given in the next lemma. 

Lemma 2.1. Let g{x,v) = g{r,Vr,vg) be a radial function on B x M^. Then g satisfies the 
specular boundary condition if and only if 



g D h dvdx = — / T) g h dvdx 

(either + or -) for all radial functions h with v-compact support that satisfy the specular condi- 
tion. 

Proof. Integrating by parts in x and v, we have 

{g D^/i + h} dvdx = 27r / gh v ■ Cr dv 

If g satisfies the specular condition, then g and h are even functions of Vr = v ■ Cr on 5ID), so that 
the last integral vanishes. Conversely, if J gh v ■ e^ dv = on (9D, then J g{l,Vr,VQ) k{vr,vg)dv = 
for all test functions that are odd in Vr, so it follows that g{l,Vr,VQ) is an even function of Vr. □ 

Therefore we define the domain of to be 

dom(D±) = [g£n\ B^gen, {B^g, h)n = -{g, B^h)n, V/i € c}, (2.10) 

where C denotes the set of radial functions h with t;-compact support that satisfy the specular 
condition. We say that a function g £ Ti with B^g € Ti satisfies the specular boundary condition 
in the weak sense if 5 € dom(D^). Clearly, dom(D^) is dense in Ti since by Lemma l2.1l it contains 
the space C of test functions, which is of course dense in Ti. 
It follows that 

{B^g,h)n = -{9,B^h)n (2.11) 

for all g,h G dom(D^). Indeed, given h € dom(D^), we just approximate it in by a sequence of 
test functions in C, and so (j2.1ip holds thanks to Lemma |2.1[ 

Furthermore, with these domains, are skew-adjoint operators on Ti. Indeed, the skew- 
symmetry has just been stated. To prove the skew-adjointness of D+, suppose that f^gGH and 
{fjh)-}i = —{g,B^h)'H for all h G dom(D^). Taking /i G C to be a test function, we see that 
B^g = / in the sense of distributions. Therefore (|2.1ip is valid for all such h, which means that 
g G dom(Z)+). 

2.4 Growing modes 

Now we can state some necessary properties of any growing mode. Recall that by definition a 
growing mode satisfies f^GH and E, B G L^(D). 
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Lemma 2.2. Let (e^V^, e^*E, e^*B) with KeA > be any growing mode. Then E,B G H^{B) and 

(l/^P + |D^/^P)^<oo. 

xI82 l/ie I 

Proof. The fields are given by (|2.2p where '(/' satisfy the elliptic system (j2.4p with Dirichlet or 
Neumann boundary conditions, expressed weakly. The densities satisfy (j2.8p . Explicitly, 

A/^ + D±/± = ±/i^i)ra^99 ± /i^i)ra^(rV') ± XifJ-t^e + (2.12) 

This equation implies that D^/^ € Ti. The specular boundary condition on is expressed weakly 
by saying that € dom(D^). Dividing by and defining = f^/\fif\, we write the equation 
in the form A^f^ + T)^g^ = h^, where the right side belongs to T-L = ±\{^ ^ 1^^) thanks to 
the decay assumption (jl.l2p on fj,^. 

Letting We = \fif\/{e+\fif\) for e > and = w^g^, we have {Xg^ + D^g^, ge)H = {weh'^,ge)H- 
It easily follows that g^ € Ti. In fact, g^ G dom(D^), which means that the specular boundary 
condition holds in the weak sense, so that (j2.1ip is valid for it. In (12. lip we take both functions to 
be ge and therefore (D'^ge, ge)'H = 0. It follows that 

|A|||c/e|||{ = \{weh^,ge)n\ < \\h^\\n\\9e\\n- 

Letting e — t- 0, we infer that (7^ G Ti, which means that ff If^l"^ /\fi^\dvdx < 00. 
Now the elliptic system for the field is 

-A^ = /(/+ - f-)dv, (A2 - A,)V = / veif' - r)dv 



together with the boundary conditions dr(p{t, 1) = 0, Vl^j 1) = 0? which are expressed weakly. 
Because of ff If^]"^ /\iJ,f\dvdx < 00, the right sides of this system are now known to be finite 
a.e. and to belong to L'^{0). So it follows that ip^^p e H'^{1D)) and E,B G H^{0). This is the 
first assertion of the lemma. Nevertheless, we emphasize that D^/^ does not satisfy the specular 
boundary condition. However, directly from ()2.12p it is now clear that ff \T)^ f^l"^ /\fi^\dvdx < 00. 
This is the last assertion. □ 

3 Linear stability 
3.1 Formal argument 

Before presenting the stability proof, let us sketch a formal proof. We consider the linearized RVM 
system (|2.8p . For sake of presentation, let us consider the case with one particle / = /"*", and thus 
drop all the superscripts +. We have the linearized equation: 

dtf + D/ = fxe^if + fMpB^rtp) + {^eve + rfXp)dtip. 

The key ingredient for stability is the fact that the functional 



dvdx + 



lEp + IBP 



dx 
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is time- invariant, which can be found by formany expanding the usual nonhnear energy-Casimir 
functional around the equilibria. Next, we then write the linearized equation in the form 

dtif - fJ-eVei' - rfipip) + D(/ - fie^p - rfipip) = 0. (3.1) 

We observe that (/ — /Xe%^ — rfipip) stays orthogonal to kerD for all time, and that in case of 
stability we would expect that / — fieV — '''fIp^p asymptotically belongs to ker D. That is, if V is the 
projection onto the kernel, one would have at large times 

/ - fie^ - rupip = V{f - /ie<^ - rfipi>), V{f - HeVe^ - rupip) = 0. 

Adding up these identities, we obtain 

/ = ^e(l - V)ip + r^pip + HeVivgilj). 

This asymptotic description of / is essential to the proof of stability, which we will prove rigorously 
in subsection 13.41 see (|3.14p . Next, by plugging the identity into the functional Z{f,ip,ip), we will 
obtain 

I{f, if, 4^) = {A^.ip, ip)L2 + iP)l2 + / \dtlP\^ dx, 

in which the operators are defined as in ()1.14p . Using the identity A^^p = B^ijj, which is obtained 
from the Poisson equation, to eliminate 99, we can formally write 

^) = i;)L2 + [ \dti^\^ dx, (3.2) 

Jo 

with := A2 + (B^)* {Ai)~^ . (In fact, we are only able to prove a reverse inequality (>), which 
however suffices for stability). Now, if we assumed that there were a growing mode of the linearized 
system, then the time-invariant functional I{f, (fi, ip) must be zero. Thus ([32]) shows that C° ^0. 
That is, > formally implies the stability. 



3.2 Key operators 

In this subsection, we shall derive the basic properties of the operators A^ and B^ defined in (jl.l4p . 
Let us recall that are the orthogonal projections of T-L onto the kernels 

ker(D±) = |/ G dom(D±) D±/ = o|, 

and the key operators: 

A21P = -Arip - j rv0{^i^ + i^ip) dvip - j vel^iJit'P^ {vei^) + f^e'P' (veip)^ dv (3.3) 
B°i> = r [ + dvil;+ [ LtV+ (ve^l^) + fi-p- {vei^)) dv. 
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We also recall that the functional space V consists of functions in ff,^(B) that satisfy the Neumann 
boundary condition and have zero average over D, whereas consists of functions in H'^^CB) that 
satisfy the Dirichlet condition on dO. 

Lemma 3.1. 

(i) Ai is self-adjoint and positive definite on L^(D) with domain V. Ai is a one-to-one map 
from V onto the set {5 G | fj^g dx = 0}. In particular, {Ai)~^ is well-defined on the range of 

(a) is a hounded operator on -L^(D). 

(ii) A^ and are self-adjoint on L^(D) with common domain . 

Proof. First, since the projections preserve the radial symmetry, the operators A^ and B'^ also 
preserve the symmetry. Next we observe that all the integral terms in (j3.3p are bounded operators 
in L^(D). For example, we have 

j^j l^ti^T^^^ dvdx <sup(^ j \^'^\dv^\\^\\L2\\V^if\\L2 < Co||V'||l2 ||v3||i2 , 

for some constant Cq that depends on the decay assumption ()1.12p on ^'^ . The other integrals 
are similar since vg is bounded by one. This proves (ii) and also proves that the integral terms 
in A^ and A^ are relatively compact with respect to —A and — A^, which have domains V and 
V^, respectively. Thus, A^ and ^2 ^'^^ well-defined operators on L^(D) with domains V and V^, 
respectively. 

Since are self-adjoint on it is clear that all three operators A^, A^, and are self-adjoint 
on L^(D). Note that the function 1 belongs to the kernel of D^. To prove the positivity of A^l, we 
use the orthogonality of and 1 — in T-L to write 

- V^)ip dvdx = - tJit\{l- V^)ip\'^ dvdx > 



since /x^ < 0. Thus Ai is a nonnegative operator, and Ai(p = if and only if 99 is a constant. 
Since A^ has discrete spectrum, it is invertible on the orthogonal set to its kernel, that is, on 
{5 € V I f^g dx = 0}. In order to prove the invertibility of Ai on the range of B^, we note that 
by the self-adjointness of in T-L, we have 

fifV^ (veip) dvdx = / / fJ't'^oip'P^i^) dvdx = / / fJ-fveip dvdx. 
JO J Jo J Jo J 

Thus 



B'^ip dx = 2J / / (r/ip + fifvg)i{; dvdx, 
j_JoJ 



which is identically zero by using the fact that d^gifJ.^] = velJ-f + rfi^. That is, B^tp has zero 
average, and so {Ai)~^ is well-defined on the range of B^. □ 
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3.3 Invariants 

First we consider the linearized energy. 

Lemma 3.2. Suppose that {f^,(p,il)) is a solution of our linearized system (j2.8p . (12. 9p . and (ll.4p 
such that G C^{R, Ll^^^^^{0 x M^)) and ip,ij G {R; H'^ (B)) . Then the linearized energy 
functional 



r 1 



dvdx + 



\dtM^ + \v^\^ + hdr{r^)\ 



dx 



is independent of time. 

Proof. For convenience of calculation, we denote the three terms as 

1 



±,1/ T'^/^pV'l -rfipVe 

|V(^|2 + |at^p + -L|a,(rv^)l 



dvdx 



dx, 



so that 



Now taking the time derivative of Zy(f'^,ip) and then using the linearized Vlasov equation ([2T 
for /+, we get 

1 d 



2dt 



1 



■(/+ - r^t^%l)){dtf^ - rii^dti^) - r^i^ve'4'dt'4' 



dvdx 



I I -^Af^ -r^l+^l,){-J}+f+ + ^itT^+ip + ^l+J}+{ri,)+^itie^tA-r^i^v^^^ 



dvdx 



rfipfi'^ipB'^ip — (/Up )^r-(/'D+(r^) 



dvdx. 



Among the nine terms, we have used the fact that two terms with dtip exactly cancel because 
/ie < 0. Some of the remaining seven terms cancel, as follows. First, we observe that the sixth 
term vanishes, upon writing D^ip = v^dr^ and noting that ^ is even in Vy. Next, by using the 
skew-symmetric property (j2.1ip of D'^, which is the weak form of the specular boundary condition, 
the first and seventh terms are also zero. Now the third and fifth terms can be combined to get 
/ip D"''(r'0/"'"), whose integral again vanishes due to the boundary condition ^ = 0. We have used 
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the fact that both ^'^ and /i^ belong to the kernel of D'^. Only the second and fourth terms survive, 
so we have 



/+D+V? - f+vedt^ 



dvdx 



'El • V dvdx, 



in which we have used the fact that 'D'^ip = Vrdr(p together with the definition E = —driver — dtipee 
Entirely similar calculations hold for We therefore obtain 

1 d 

2diK 



^^(x+(/+,V)+X^(r,V')) = / I E.v{f+-r)dvdx = I^B.jdx. (3.4) 

dx. Taking its time derivative and 



Next, by ()2.ip we may write Xa/((/9, -0) = /j 
using the Maxwell equations, we obtain 



Er + B 



B-dtE + B- dtB 



dx 



dx. 



- E • j + E • (V X B) - B • (V X E; 
/ E • j + / (E X B) • n(x) dS^ = - / E • j 

Jd JdD Jo 



dx 



Here we have used the fact that (E x B) -n = (Ex n)-B, which vanishes due to the perfect conductor 
boundary condition E x n = 0. Together with (j3.4p . this yields invariance of I{f^,ip,ip). □ 

Furthermore, we also obtain the following. 

Lemma 3.3. Suppose that {f^,(p,tj;) is a solution of our linearized system (j2.8p . (j2.9p . and ()1.4p 
such that /± € C^(M, (D x R"^)) and ip,Tp £ {R; H"^ (B)) . Then the functional 

^^(/^,V')=/ j (f^^iitvB^Tr^Agdvdx (3.5) 



are independent in time, for all g G kerD^ and for both + and —. In particular, for g = 1 in (|3.5p . 
the integrals /p /]g2 f^{t,x,v) dvdx are time-invariant; that is, the total masses are conserved. 

Proof. Parenthetically, we remark that such invariants }C^{f^, ip) can easily be discovered by writ- 
ing the Vlasov equations as in (j3.ip . Indeed, writing the Vlasov equations in the form p.ip and 
using the skew-symmetry property of as in (|2.1ip . we have 

d 



dt{f^T fJ-t^etp =F rtJ,p'ip]g dvdx 



— f^ ± fifip i: rfi^ipjg dvdx 

- ± fif(p ± rfi^ip j D^g dvdx = 



due to the specular conditions on and g and the evenness in Vr of /U^. Now, if we take g = 1 
in ()3.5p and note that dvgfi^ = fJ'tve + rfi^, the integrals fjj^ f^2 f^{t,x,v) dvdx are therefore 
time-invariant. □ 
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3.4 Minimization 

In this subsection we prove an identity that will be fundamental to the proof of stability. It involves 
the functional 

Mf^J~)-= [ [ T^lf-^ -rfi;i;\Uvdx+ [ [ -^\r + rf,;i;\Uvdx+ [ 
Jo J \^J■e\ JdJ lAiel Jn 

where if = (p{r) € V satisfies the Poisson equation 

-Aip= /(/+ - /-) dv, [ ^dx = 0, dMl) = 0. (3.6) 
J Jo 

For each fixed ip € L^(]D)), let J-^, be the space consisting of all pairs of measurable functions 
{f^,f~) depending on {r,v) which satisfy the constraints 

^\ \^ dvdx < +00, (3.7) 



and 

(^/+ - - ?^/^p V') dvdx = 0, j j (^f~ + n've^j + rfi~il?j g~ dvdx = 0, (3.8) 

for all € kerD^. Similarly, let be the space of pairs {f~^,f~) satisfying (13. 7p and 

f'^g-^dvdx = 0, / f~g~dvdx = 0, Vy^^ekerD^. (3.9) 



Note that the constraints in ()3.8p with g = I imply that for such a pair of functions there is 
a unique solution ip £ V the Poisson problem ()3.6p . In particular, is radially symmetric since 
are radially symmetric. Thus the functional J'^ is well-defined and nonnegative on T^, and its 
infimum over J-"^ is finite. We next show that it indeed admits a minimizer on J^^. 

Lemma 3.4. For each fixed ip € L^(B), there exists a pair of functions f^ that minimizes the 
functional J'^ on T^. Furthermore, 

Mf*J*) = ((e°)*(^;)-isV,V)L2 - / / /.+ |P+(i),^)|2 dvdx - f I iJi-\V-{ven^ dvdx. 

Jo J Jo J 

(3.10) 

Proof. Take a minimizing sequence in F^ such that Jip{fn ■, fn) converges to the infimum of J^. 
Since {fn} are bounded sequences in L^^|^_|_|, the weighted space associated with the constraint 

()3.7p . there are subsequences with weak limits in L^^|^±|, which we denote by f^. It is clear that 

satisfy the constraints (|3.8p . and so they belong to F^. That is, (/+, /~) must be a minimizer. 

In order to derive the identity (|3.1Up . let the pair (/*,/") G be the minimizer and let 
V?* G H'^.ijii) be the associated solution of problem (j3.6p with f^ = f^. For each {f^,f~) € F^, 
we denote 

h'^ := f+ - H^vgip - rfi+7p, h~ := f^ + n^veip + rf^ipijj. (3.11) 
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In particular, hf := ^'^VQil^^rfi^il). It is clear that (/"'", / ) S if and only if h ) G Tq. 
Since d^g [/i^] = ^i^vg + r/i^ , we have 

(/+ - /-) dv = j {h+ - h-) dv. 

Thus if we let be the solution of the Poisson equation (|3.6p . is independent of the change of 
variables in (|3.1ip . Consequently, (/^, /") is a minimizer of on J^^ if and only if {hf, h~) is a 
minimizer of the functional 



Jo{h+,h- 



j f , \h'^ + fi'^vgipl'^ dvdx + f f -——\h — fi^ vgipl'^ dvdx + f |V<y9p dx, 

Jo J IfJ^e Jo J IfJ-e Jo 



IfJ'C I Jo J IfJ-e 

on Tq. By minimization, the first variation is 

/ I -^{ht + l4veil^)h'^ dvdx + j j -^{h; - ii~ve'4>)h^ dvdx + I Vip^ -Vip dx = {), (3.12) 

Jo J fJ-e Jo J /Ue Jo 

for all {h^ , h~) € J-Q where (p solves (|3.6|) . By the Neumann boundary condition on ip, we have 

/ • Vip dx = — ip^:/S.Lp dx = / (p^^h^ — h~) dvdx. 
Jo Jo Jo J 

Adding this to the identity (|3.12p . we obtain 

/ I —rih^ + lJ.eV9ip — iJ.'^'p,^) dvdx + [ [ ^{h~ — fi'vgilj + ij,~ip^,) h~ dvdx = (3.13) 
JoJ fJ'I JoJ Aie 

for all {h^ , h~) (z J-q. In particular, we can take /i~ = in (I3.13|) to get 

-1 



{ht + — fJ-t^*) dvdx = 

Ate 

for all G ^\/\^+\ satisfying Jp J h'^g'^ dvdx = 0, for all g(+ G ker D"*". 

We claim that this identity implies hf + fJ-t'^'e'^P — fJ-t^* € kerD'''. Indeed, let 

Using the inner product in = _|. , we have 

{k^J)^ = V£G(kerD+)^. 

Because D"*" (with the specular condition) is a skew-adjoint operator on 7i, we have k^, G (ker D"*") 
kerD+. Thus 



±± 
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This proves the claim. Similarly D {/^ + rfip i/j + fj,^ ip^} = 0. Equivalently, 

On the other hand, the constraints p.Sp can be written as V^[f^ — fJ-tveip — rfi^tp) = and 
V~{f~ + /i~?}6)V + V') = 0. Combining these identities, we have 

_ _ _ _ _ _ (3-14) 

Thus, using the orthogonality of and (1 — V^) in H, we compute 

^\ft TTfi^tpl'^ dvdx = - [ [ fif\{l-V^)^^\'^ dvdx- [ f i4\'P^{v0^)? dvdx 
fJ-e JnJ Jo J 

— V^)ip^ dvdx — I I ^^|P^('06i^)P dvdx. 



Inserting these identities into the definition of J^{ft, f* )i together with the fact from the boundary 
conditions that Jp |V(/9*p dx = — ip^Atp^ dx, we obtain 



J^ift, /* ) = / - ^f* - / (1 - 'P^)^^* dv - I fi^il-V )ip^ dv 
Jo ^ J J 

fitlV^ivei^)]"^ dvdx - / / fi-lViveij)]'^ dvdx. 



dx 



By the definition ()3.3p of Ai, the first group of integrals simply equals {Ai(p>t, 

Thus it remains to prove that ^iV?* = B^ip, because Ai is invertible on the range of so 
that 93* = {Ai)~^B^il). Indeed, we plug the identities (|3.14p into the Poisson equation (|3.6p for 99^., 
resulting in the equation 

-Aip^ = / /i+(l - dv+ / - V~)ip* dv 

+ r {fi^ + fip) dvtp + / i^fitV^ivei^) + P {vei')j dv. 

In view of the definitions of Ai and B^, this identity ()3.15p is equivalent to ^Jv'* = B^^p, as 
desired. □ 



3.5 Growing modes are pure 

In this subsection, we show that if (e'^*/^, e'^*E 
then A must be real. See subsection 12.41 for the 
splitting method in [T5] to show that A is real. 



, e^^B) with SReA > is a complex growing mode, 
properties of a growing mode. We now follow the 
Let /ev and /od be the even and odd parts of / 
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with respect to the variable Vr- That is, we have the spHtting: / = /ev + /od) and furthermore, by 
inspection from the definition (j2.7p . the operators map even to odd functions and vice versa. 
We therefore obtain, from the Vlasov equations (|2.8|) . the spht equations 

f A/+ + D+/+ = \{^ltve + r/x+)V = AO,, 

By Lemma 12.21 we know that Jj \f'^\'^ /\^Lf.\dvdx < oo. It follows that the same is true for f^, so 
from the first split equation we also have j ^ \D~^ f^^\'^ /\iif.\dvdx < oo. The split equations imply 
that 

(A^ - D+')/+ = X^itB+^ - A/i+D+C^eV'). (3.16) 

Let be the complex conjugate of /+. By the specular boundary condition on in its weak form 
()2.1ip . it follows that f^^ satisfies the specular condition. (Formally, f^^ vanishes on the boundary 
9D.) However, since D~^f^^ is even in the variable Vr, f^^ also satisfies the specular condition. 

Thus when we multiply equation (j3.16p by /^d/l/^eM and integrate the result over B x M^, we may 
apply the skew-symmetry property (j2.1ip of D"*". We obtain 

A'/ [ T^\fofdvdx+ [ [ J^\B+f+fdvdx = - [ [ (xB+^Jt^ + XveiPD+JtAdvdx. 
Similarly for /~ we obtain 

A' / / -^\f-d\^ dvdx + [ f ^\B-f-^\Uvdx= f f (xB-^J-^ + Xve^^^B-J~^ dvdx. 

Jo J IfJ-e I JD J l/ie I JoJ^ ' 

Adding up these identities and examining the imaginary part of the resulting identity, we get 

25ReASmA / / (^|/+ p + 1 |y^-j2\ ^^^^ 

VI;,, I I;,, I ^ } ^ 

^^'^^foA - D~93/od) 'J'vdx - 9m / / Af;0V'(D^/od - D"/od) dvdx. 

Jo J 

Let us now use the Maxwell equations to compute the terms on the right side of ()3.17p . First, 
we recall that the second equation in ()2.4p is 

Xdr^ = I {),(/+ - /-) dv = J VriC - /o"d) dv. 

Together with the definition of D^, this yields 

-X [ [ (D+(^7+ - B~^J~^) dvdx = -X [ dr^( [ vriltd - lod) dv)dx = -|A|2 / \dM''dx, 
Jo J Jo ^ J ^ Jo 

whose imaginary part is identically zero. 
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Secondly, using the Vlasov equation for f^, we estimate 
-X [ [ ve^{B-^7td - D"7od) dvdx 



By (j2.4p the first term on the right equals 



Jo 



-Ar + X")Hx=\X\'^ I ^Idriri^)]"^ dx + X^\X\'^ / |Vp dx. 



Here we integrated by parts and used the Dirichlet boundary condition on ip. 
Putting these estimates back into (I3.17p . we obtain 

25ieA9mA / / (-^\f+^\'^ + -^\f^^\'^)dvdx = -2^eXQmX\X\'^ [ IV'P dx. 

Jn J ^ l/Ue I l/Ue I ^ Jo 

The opposite signs of the integrals imply that A must be real. 
3.6 Proof of stability 

With the above preparations, we are ready to prove the following stability result, which is half of 
Theorem II. li 



Lemma 3.5. If > 0, then there exists no growing mode (e"***/^, e'^'^E, e'^*!?), with 5ReA > 0. 

Proof. Assume that there were a growing mode (e^*/^, e'^*E, e^*i?). For the basic properties of 
any growing mode, see Lemma 12.21 By the result in the previous subsection, it is a purely growing 
mode, and thus we can assume that (/^, E, B) are real- valued functions. By the time-invariance in 
Lemma [3.2l the functional ip, ip) must be identically equal to zero, where tp and ^ are defined 

as usual through the relations ()2.ip . That is. 



I{f^,ip,^P) = Mf^,f-)- I I rvei^ + iJi^m^ dvdx + I \X\'\il>\' + ^\dr{r 



|2l;.|2 I 1 lo t^j.\\2 



dx = 0. 



Furthermore, all the expressions lC^{f^ ,4)) defined in (13. 5p must be zero. The vanishing of the 
latter integrals is equivalent to the constraints in (j3.8p and therefore the pair {f~^,f~) belongs to 
the function space J^^. We then apply the Lemma [33] to assert that Jxi,{f~^,f~) > Ji^ift ■• f*)- 
Thus we have 

X(/±,</.,V') > ((S°)*(^?)-^eV,V')L2 - / / ^it\V^{ven^ dvdx - [ [ f,~\V-{ven^ dvdx 

dx. 



rveifi; + fi^m^ dvdx + 1^ |A|2|V;|2 + l|9,,(rV')|' 



(3.18) 
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In addition, from (|3.3p an integration by parts together with the Dirichlet boundary condition on 
ip yields 



fJ't\'P^M)\'^ dvdx - fj.~\V~{vgip)\'^ dvdx. 
Jd J 

Putting this calculation into (j3.18p and using the definition of C^, we then get 



This is obviously a contradiction since we are assuming > 0. Thus there exists no growing mode 
for the linearized system. □ 

4 Linear instability 

We now turn to the instability part of Theorem 11.11 It is based on a spectral analysis of the 
relevant operators. We plug the simple form (e'*'*/^i c'^'^E, e^*i?), with A > 0, into the linearized 
RVM system (j2.8p to obtain the Vlasov equations 



(A + D+)/+ = /z+D+^ + /i+D+(rV^) + A(Ai+*e + rAi;)^ 
(A + D^)/" = -M;rD~(^-^~D~(rVj) - A(^;:ue + r/i")V' 



(4.1) 



and the Maxwell equations 

-A^ = j{f+-r)dv, 



(4.2) 



(-A, + A2)V' = j ve{f^-f-)dv. 

As before, we impose the specular boundary condition on the Neumann boundary condition 
on if, and the Dirichlet boundary condition on ijj. Recall that {f^,Lp,il)) is a perturbation of the 
equilibrium. 

4.1 Particle trajectories 

We begin with the + case (ions). For each {x,v) S D x M^, we introduce the particle trajectories 
{X~^ {s;x,v),V^ {s]x,v)) defined by 

X+ = T/+, y+ = E°(X+) + T/+ X B°(X+), (4.3) 

with initial values (X"'"(0; u), 1/^(0; x, f)) = {x,v). Because of the regularity of E*^ and B*^ in 
D, each trajectory can be continued for at least a certain fixed time. Thus each particle trajec- 
tory exists and preserves cylindrical symmetry up to the first point where it meets the boundary. 
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The trajectories reflect at the boundary 9B at most a countable number of times as s — > ±00. 
Let So be any point at which the trajectory X^{so—;x,v) belongs to dlDi. In general, we write 
/i(sib) to mean the limit from the right (left). By the specular boundary condition, the trajectory 
{X~^{s]x,v),V^{s;x,v)) can be continued by the rule 

(X+(so+; X, v),V-^{so+; X, v)) = (X+(so-; x, v), V+{so-; x, v)), (4.4) 

with V = { — Vr,Vg). Thus X'^ is continuous and has a jump at sq. Whenever the trajectory 
meets the boundary, it is reflected in the same way and then continued via the ODE (|4.3p . Such 
a continuation is guaranteed for some short time si past sq by the standard ODE theory. Since 
EO and are smooth in D, the additional time |si — sq\ is bounded below by some fixed 
positive time Tq independent of sq. This shows that the trajectories can bounce at the boundary 
at most a countable number of times as \s\ — > 00. When there is no possible confusion, we will 
simply write {X~^{s), V~^{s)) or {R'^{s), Vj^{s), Vf^{s)) for the particle trajectories. The trajectories 
(X" {s),V~ (s)) for the — case (electrons) are defined similarly. 

Lemma 4.1. For each {x,v) G B x M^, the particle trajectories {X^{s;x,v),V^{s;x,v)) are piece- 
wise smooth m s € M, and for each s G R, the map {x,v) i-7> {X^{s;x,v),V^{s;x,v)) is 
one-to-one and differentiable with Jacobian equal to one at all points {x, v) such that x dO and 
X^{s;x,v) 9B. In addition, the standard change-of-variables formula 

/ / gix,v)dxdv= / g{X^{-s;y,w),V^{-s;y,w)) dwdy (4.5) 

is valid for each s € M and for measurable functions g for which the integrals are finite. 

Proof. For each {x,v), the particle trajectory {X^{s;x,v),V^{s;x,v)) is smooth except when it 
hits the boundary dD, which happens countably many times. So the first assertion is clear. Given s, 
let S be the set {x,v) in B x such that X^ {s;x,v) 9B. Clearly, S is open and its complement in 
B X has Lebesgue measure zero. For each s, the trajectory map is one-to-one on S since the ODE 
()4.3|) and (14.40 are time-reversible and well-defined. In addition, a direct calculation shows that the 
Jacobian determinant is time-independent and is therefore equal to one. The change-of-variable 
formula (j4.5p holds on the open set S and so on B x M^, as claimed. □ 

Lemma 4.2. Let g{x, v) be a radial function on B x M^. If g is specular on 5B, then for all s, 
g{X^{s;x,v),V^{s;x,v)) is continuous and also specular on 9B. That is, 

g{X^{s;x,v),V^{s;x,v)) = g{X^{s;x,v),V^{s;x,v)), V {x,v) e dB x M^, 

where v = {—Vr,VQ) for all v = {vr,vg). 

Proof. It follows directly by definition ()4.3p and ()4.4p that for all x € 9B, the trajectory is unaffected 
by whether we start with v or v. So for all s we have 

X^{s;x,v) = X^{s;x,v), Vg^{s;x,v) = Vg^{s;x,v), \V,^{s;x,v)\ = \V,^{s;x,v)\. (4.6) 

In fact we have V^^{s; x, v) = Vj^{s; x, v) for any s at which X^{s; x, v) (9B, while V^{s-\-; x, v) = 
—V^{s—;x, v) for s at which X^{s; x, v) G (9B. Because g is specular on the boundary, it takes the 
same value at Vr and —Vr- Therefore g{X^[s),V^{s)) is a continuous function of s at the points 
of reflection. It is specular because of the rule 1)4.40 . □ 
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4.2 Representation of the particle densities 

We can now invert the operator (A + D^) to obtain an integral representation of By definition 
of the operator from (j2.7p and the trajectories (X+(s), V~^{s)) from (j4.3p and (j4.4p . we have 

° e''B^g{X+{s),V+{s))ds= r e''^g{X+{s),V+{s)) ds 

= g{x)- r \e^^g{X+{s),V+{s))ds, 



for functions g = g{x,v) which belong to dom(D'^). Multiplying the Vlasov equations (14. Ih by 
c*"* and then integrating along the particle trajectories from s = — oo to zero, we readily obtain 

rO rO 

f+{x,v) = ntv^ + rfJ-^i^- fJ-t Xe^'^{R+{s)) ds+ Ae^'/x+y+(s)'i/'(i?+(s)) ds. 

J —oo J —oo 

A similar derivation holds for the — case. For convenience we denote 

rO 

As 

At 

oo 



Qf{g){x, v) := / Xe^'g{X^{s; x, v), V^{s; x, v)) ds. 



In particular, by Lemma 1^21 Q^ig) is specular on dO if g is. Thus we have derived the integral 
representation for the particle densities: 

v) = fitil- Qx)v^ + rfi^i; + fitQxivei^) (4.7) 

4.3 Operators 

We next substitute (|4.7p into the Maxwell equations (|4.2p . We introduce the operators 
A^ip : = -Aip - J - Q+)v9 dv- J - Q'^)ip dv, 

r r (^-^^ 

B^^: = - Ut{l- Qt){v0i^) dv - / ^-(1 - QDiveij) dv, 



(B^r^P : = - j - dv- j - Q~x)^ dv. 

We also introduce 

■= + {B^y{A\y^B^. (4.9) 

We then have 

Lemma 4.3. The Maxwell equations (j4.2p are equivalent to the equations 

A\ip = B^i}, + (-S^) V = 0. (4.10) 
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Proof. Using (j4.7p . we write the Poisson equation for ip as 
-Aip = j {f+-f-){r,v) dv 

= j - Qt)f dv + j /i^(l - ^'"^ j ^^^P )^^^ 

+ j fitQtM) dv + j fi-Q-{vg^) dv. 

Note that d^g [fi'^] = rfi^ + vgfif so that 

j riii'^ + ii~)dv = - J veifit + l^e) dv. 
This gives the first equation in (j4.10p by definition. Similarly we write 
{-Ar + X^)^P = j veif^ - r)ir,v) dv 

= j vefi'^{l - Ql)ip dv + j ven~{l - Q'^)ip dv + j rved^ip + fi~)dv'ip 
+ / vgfi^Q'livgij) dv + / vel^eQ-xi^eip) dv, 



which is equivalent to the second equation in (j4.10p . □ 

As in Lemma |3. 11 we now state some properties of these operators. We recall that the spaces 

V and V^, which are defined in Section [1.21 incorporate the boundary conditions. 

Lemma 4.4. For any A > 0, 

(i) Ai is self-adjoint and positive definite on Li^-iji) with domain V. Moreover, Ai maps from 

V one-to-one onto the set 1^ := {g ^ : j^^g dx = 0}. 

(ii) is a hounded operator on L^(D) with its adjoint operator {B^)* defined in (|4.8p . The 
range of B^ is contained in 1-*-. 

(Hi) A2 and are self-adjoint on L^(D) with their common domain V^. 

Proof. We first check the self-adjointness of A'j and the formula for {B^)* . Since ji^ is constant on 
trajectories and in view of (j4.8p . it clearly suffices to prove, for smooth functions g and h specular 
on the boundary, that 

,tH.,.)QH,i,,.)) = / I ,t9(^.i,)QtiM.M ivi., (4.11) 

Jn J 

where we denote v = {—Vr,vg). In order to prove (|4.1ip . we recall the definition of and use the 
change of variables 

iy,w) := {X+{s;x,v),V+{s;x,v)), {x,v) := iX+{-s;y,w),V+{-s;y,w)), 
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which has Jacobian one where it is defined (see Lemma l4.ip . So we can write the left side of (j4.1ip 



as 

rO 



\e^^ lX^h{x,v) giX'^ {s;x,v),V^ {s;x,v)) dvdxds 

■0 



Ae ''^lJh{X+{-s]y,w),V+{-s■,y,w)) g{y,w) dwdyds. 

Observe that the characteristics defined by ()4.3p and ()4.4p are invariant under the time-reversal 
transformation s ^ —s, r r, Vr ^ —Vr, and vg ^ vg. Thus 

X^{—s; X, v) = X~^{s; x, v), V'^{—s; x, v) = V~^{s; x, v), 

at least if we avoid the boundary. Using this invariance, we obtain 

Xe^"^ fj,'^ h{x , v) g{X~^{s;x,v),V~^{s;x,v)) dvdxds 
rO r r 

Xe^''n^h{X+{s;y,w),V^{s;y,w)) g{y,w) dwdyds 
Xe''^'^fi'^h{X~^{s;x,v),V~^{s;x,v)) g{x,v) dvdxds, 



-oo 

° ' ' A, 

L I ^'^ 



in which the last identity comes from the change of notation {x,v) := {y,w). By definition of Q^, 
this result is precisely the identity (I4.1ip . A similar calculation holds for the — case. This proves 
the adjoint properties claimed in the lemma. 

Next we show that all the integral terms in (|4.8p are bounded operators on L^(D). For instance, 
we have 

fiti'Qt'P dvdx = I / Xe^'^fi^^pipiX+is)) dvdxds 
< ( f Xe^' [ [ IfitM^ dvdxdsV^^ f Xe^' [ [ \fitMX-^{s))\^ dvdxdsV^^ 
<sup(| \iit\dv)m^.J^\\j^.^. 

(4.12) 

In the last step we made the change of variables {x,v') = {X^{s;x,v),V~^{s;x,v)) in the integral 
for (f, which is possible thanks to (I4.5p . Similar estimates hold for the other integrals since vg is 
bounded by one. This proves that is bounded on L'f.[n) and also that the integral terms in Ai 
and A2 are relatively compact with respect to —A and — A^, respectively. Therefore and A2 
are well-defined operators on L^(D) with domains V and V^, which are the same as the domains of 
—A and — A^, respectively. 

Taking ■0 = in the previous estimate, we have 



fJ'tv^Qxf dvdx 



00 



< / Ae^"* / / \h^\\lp\^ dvdx = - / / n'^\ip\^ dvdx 
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so that 

- / / ^it^ (1 - Ql)^ dv > 0. 
Jo J 

Thus Ai > 0, and Ai^p = if and only if (/9 is a constant. Since Ai has discrete spectrum, it is 
invertible on the set orthogonal to the kernel of Ai- That is, it is invertible on G V : jj^g dx = 
0}. For the invertibility of A\ on the range of B^, we note by ([ilT]) and Q^(l) = 1 that 

I /"/^^(l - Qt){vei^) dvdx = [ [ fifve^l - Qx)il) dvdx = 0. 
Jo J Jo J 

This shows that B^ip G {<? € V : fj^g dx = 0} for all ■0. Thus {Ai)~^ is well-defined on the range 
of B^, and so is well-defined. The self-adjoint property of C"^ is clear from that of Ai- □ 

Part (i) of Lemma 14.41 in particular shows that for each ijj € L^(D) there exists a unique radial 
function € H^{0) that solves 



Aif = B^ip, / (pdx = 0, dr<f{l) = 0. (4.13) 
Jo 

4.4 Construction of a growing mode 

Lemma 4.5. // ^ 0, i/ien i/iere exisi a X > and a nonzero function tp G H'^'^{D) such that 
C^ijj = and satisfies the Dirichlet condition on the boundary. 

Proof. The proof is similar to the one and a half dimensional case given in |15j so that we merely 
outline the main steps as follows. 
(i) > for large A. 

(a) For all ip G L"^, C^ip converges strongly to C^ijj in as A ^ 0, and thus ^ ^ when A is 
small. 

(Hi) The smallest eigenvalue := inf^^ {^^"ip, of is continuous in A > 0, where the infimum 
is taken over iPeV^ with \\iP\\l2 = I. 

These three steps imply that k'^ must be zero for some A > 0, from which the lemma follows. 
To prove (i), it is easy to see that is nonnegative for large A, since (B^)* {Ai)~^ B'^ > and 
{A2'ip,ip) 1,2 is sufficiently large when A is large. As for (ii), to show the convergence of to as 
A ^ 0, we use the remarkable fact, proved in [HI Lemma 2.6], that for all 5 G the strong limit 

^Inn Q^g = V^g (4.14) 

is valid in the +, = 1-L norm. Here the are the orthogonal proiections of +, onto the 

kernels of D^. However, it should be noted that the convergence is not true in the operator norm. 
For all G L^, we use (|4.14p and write 



dv, 
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thereby obtaining the convergence of .4,2^ to A21P in L^. Similarly, Ai and B''^ also converge 
strongly in to Ai and B^, respectively, and so does to C^. Finally, estimates similar to (j4.12p 
yield 



< 



We 



Xs 



fie 



I 



dvdx 



1/2 



\fiJ\\i;{X+{s))\^ dvdx) ds 



1/2 



<Co 



lAe 



Xs 



fien ds) \ml, < Col log A - log/xlllV'lli^, 



and thus 



(>l^^-^^V,V'> < Co(\X-fi\ + |logA-log/i| 



i2 



for all A,/i > and V € -L^. Similarly, we obtain the same estimate for C"^, which proves the 
continuity of the lowest eigenvalue k"^ of >C^. □ 

Using ip, we can now construct the growing mode. 

Lemma 4.6. Let X^ip he as in Lemma \4-5l let tp be as in (|4.13p . and let be defined by (I4.7p . 
Then (e^*/^, e^V, e^*^) 

is a growing mode of the linearized Vlasov- Maxwell system. 

Proof. Because C^il) = and due to the definition of 93, both parts of (j4.10p are satisfied. Therefore 
(j4.2p is satisfied. These are the first and third Maxwell equations in (j2.4p together with the boundary 
conditions for ip and ip. Next, the specular boundary condition for follows directly by definition 
()4.7p and the fact that Q^{g) is specular on the boundary if g is. It remains to check the Vlasov 
equations and the middle Maxwell equation in ()2.4p . namely Xdr(p = jr- 

We begin with the equation for /"•". Recall that {X~^{t; x, v), V^{t; x, v)) is the particle trajectory 
initiating from {x,v). Evaluating (j4.7p along the trajectory, we have 

f+{x+{t),v+{t)) = + - /.+ r Xe'Mx+i^;X+{t),v+m ds 



Ae^^y+(s;X+(t),y+(t))V(X+(.;X+(t),y+(t))) ds. 



By the group property {X+{s;X+{t),V+{t)) = X+{s + t) and V+{s;X+{t),V+{t)) = V+{s + t), 
together with integration by parts in s, we have for each t 



f+{X+{t),V+{t)) = fit^{X+{t))+4R+{tmx+{t)) 

fO 



Xe 



Xs 



99(X+(s + t)) - V+{s + t)il;{X+{s + t)) 



ds 



/.+i2+(t)V(X+(t)) + M+e 



-At 



dMX+{s)) + XV+{s)^P{X+is)) 



ds. 
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Differentiation of this identity yields 

^^[e''f+{X^{t),V+{t))) = ^,;^^(^e''R^{t)i;{X+{t))) + pte''[^^^^ 

We evaluate the above identity for t G (0, e) and let e ^ 0. Note that by Lemma l4.ll the functions 
f+{X+{t),V+{t)), ip{X+{t)) and ^p{X+{t)) are piecewise smooth. Using the evolution K3\\ 
and (|4.4|) . we obtain 

This is the Vlasov equation (j2.8p for f~^. A similar verification can be done for /~. 

Finally, we verify the remaining Maxwell equation Xdr^p = jr- Indeed, by performing the 
integration in v of the Vlasov equations (|2.6|) . we easily obtain + V • j = 0. Together with the 
Poisson equation in (j4.2p . this yields 

- (^dr + (Xdrf) = -XAif = Xp= -(j3r + ^) jr 

Thus r{Xdr(p — jr) must be a constant. However, at the boundary r = 1 we have dr(p = and 
jr = by the specular boundary condition on So Xdr^p — jr = 0. □ 



This completes the proof of Theorem 11.11 

5 Examples 

The purpose of this section is to exhibit some explicit examples of stable and unstable equilibria, 
and thereby prove Theorem 11.21 

5.1 Stable examples 

By Theorem 11.11 the condition for spectral stability is 

£° = ^^ + (S°)*(^?)-ie° > 0. (5.1) 

For each ip in the domain of (thus in particular satisfying the Dirichlet boundary condition), we 
have 

where ip solves Aiip = B^ip with the Neumann boundary condition on (p. Recall that 
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in which jj^ denote fi^{e^,p^) = fi'^{{v) ± ip^,r{vg Integration by parts, together with the 

boundary conditions for (p and ip, and the orthogonahty of and 1 — (separately for + and 
— ) lead to the expressions 



{A'i^,ip)L2= / idM^dx- / fit\{i-r^){v)\^ dvdx- / U;\{i-v-){ip)\^ dvdx, 

Jb Jo J Jo J 

=^(|a,^|2 + l|V^|2) dx- rve{4 + ^i-p) dv)\^\^ dx M 



dvdx. 



Due to the assumption /z^ < 0, it is clear that (^iV^, <^)l2 > and so are the first and last terms 
in (^2^,^^)1,2. We now exhibit two explicit sufficient conditions for {j^ijj^ijj) ip. to be nonnegative. 
This is Theorem 11.21 (i) and (ii). 

Theorem 5.1. Let {ii^,ip^,ip^) be an inhomogenous equilibrium. 
(^)If 

Pl^p{e,p)<0, y e,p, (5.3) 

then the equilibrium is spectrally stable provided that ip^ € and ip^ is sufficiently small in L°° . 
(^^) If 

\4{e,p)\ < Y^^, (5-4) 

for some 7 > 2, with e sufficiently small and ip^ = but ijj^ not necessarily small, then the 
equilibrium is spectrally stable. 

Proof. First consider case (i). We only need to show that > 0. Let us look at the second 
integral of {A2ip,'ip)]^2 in (15. 2p . By the definition p.9p of p^, we may write 



dx. 



rvQfj.p{e ,p ) dv = J (v) p^fj.p{e ,p )dv^rip'J{v) fj,p{e ,p ) dv, 
in which the first term on the right is nonpositive due to ()5.3p . Therefore we have 
~ j {j ^^silJ-p + IJ-p) dv^ dx> j j (■i;)~^(/j,+ - dv rV'^lV'P dx 

> -sup|^0|(sup f {vrWii;\ + \ii~\)dv) I r\ 

r ^ r J ''Jo 

Now by the Poincare inequality, 

/ riV'P dx<co [ + ^IV'P) dx, 

Jo Jo^ J 

for some constant cq. In addition, thanks to assumption (jl.l2p . the supremum over r € [0,1] of 
/ I + I) dv is finite if 'pP is bounded. Thus if the sup norm of ■0'^ is sufficiently small, 
or more precisely if V'^ satisfies 

cosup|0°|(sup J{vyWfi+\ + \fi-\) dv^ <1, (5.5) 
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then the second term in (^2^') '0)l2 is smaller than the first, and so the operator A2 is nonnegative. 

Case (ii) is even easier. As above, we only have to bound the second term in {A2'4' , i^) ■ Using 
\rvg\ < 1 and e = {v) together with (|5.4|) . we have 

If e is sufficiently small, the second term is smaller than the positive terms. □ 
5.2 Unstable examples 

For instability, it suffices to find a single function in the domain of such that (C^ipjip) ^2 < 0. 
We shall construct some examples where this is the case. We limit ourselves to a purely magnetic 
equilibrium (/x=^,E°,S°) with = and = i9^(r'0°). Thus e = {v) and = r{ve ± 1/^°). In 
this subsection, we shall also make the assumption that 

lj,+ {e,p) = fi-{e,-p), Ve,p. (5.6) 

This assumption holds for example if /x"*^ = /i~ is an even function of p. It greatly simplifies the 
verification of the spectral condition on C^. 

We now show that assumption (|5.6|) implies that the operator vanishes and so simply 
reduces to Indeed, let us recall that 

B^i; = r I (;U+(e,p+)+/x-(e,p-)) dv^l^ + ^ (^f^t {e,p+)V^ (veij) + f^- {e,p-)V- {vei^)^ dv, 

in which e = {v) and p^ = r{vQ it ip^). For the first term in B^tp, we again note by (|5.6p that the 
function 

A^p (e,P"^) + (e,p") = -fip{{v),-r{ve + + fip{{v),r{ve - -0°)) 
is odd in vg. Thus, the first integral in B^ vanishes. As for the second integral, we note that 

D- = Vrdr + + ^i-ve){-ve))d,^ - (i),^^ + ^Vri-vg))d-,,,. 

That is, D~ acting on functions f{vr,—vg) is the same as acting on f{vr,V0) (ignoring the 
dependence on x). As a consequence we have 

V^{f{Vr,Ve)){Vr,Ve) = V~{f{Vr, -Ve)){Vr, -Vg). 

Using this identity together with the fact that V^{—f) = —'P^{f), we have 

fi^{e,r{vg + il/^)!^^ {v0'4)){ve) + fj.~{e,r{ve - il)^)!^' {veip){ve) 
= - fi-{e,r{-ve - iP^))V- {vei^){-ve) - fit{e,r{-ve + iP^))V+(vei^)(-ve) 
= -fit{e,ri-ve + 4^°))V+{vetp){-ve)-fi-{e,r{-ve-ip^))V~(vei^)(-ve). 
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Thus the function fi'^V~^{voip)+fi^ V (vetp) is odd in vq, so that the second integral in B^Tp vanishes. 
Similarly, we easily obtain 

A21P = -Arip - 2 j rv0fj.p{e,p^) dvip -2 j ve^j," {e,p~)V~ {vei^) dv. (5.7) 

We summarize the above considerations in the following lemma. 
Lemma 5.2. If {ij,^,0,ip^) be an equilibrium satisfying (|5.6p . then = and 

= Al 

for ^2 (j5.7p . In particular, A2 ^ implies the spectral instability of {fi^,0,tp^). 



5.2.1 Homogeneous equilibria 

We start with the homogenous case EO = and B° = 0, in which case the linearized Vlasov 
operator reduces to D = • = Vrdr- The projection V = is simply the average 

V{ip) = — / 'ip{r) dx = 2 / ijj{r)rdr 
71" Jo 

for any radial function ip = ip{r). In addition, noting that T){rvg) = 0, we have 



V{veip) = rvgV{—) = rvg- I r ^^p{r) dx = 2rvgiljR, tpR 



vr 



^(r) dr. 



Thus, as in ()5.2p . we obtain the basic identity 

=^(|9.^|' + ^IV'P) dx-2 
+ dx-2 




rvg^p dvj IV'P dx — 2 
{v)^^pfip dvj IV'P dx - 



/Ug \'P{vgip)\'^ dvdx 
r'^V0i-i~ dvdx\ipji\'^ 

= 1 + 11 + III, 
in which p = rvg. 

Theorem 5.3. Let fi^ = fi^{e,p) be an homogenous equilibrium satisfying (j5.6p . Assume also that 



PfJ-p (e,p) > cop^iy{e), 



V e,p, 



(5.8) 



for some positive constant cq and some nonnegative continuous function i/(e) such that 1/ ^ 0. 
Then there exists a positive number Kq such that both of the resettled homogenous equilibria (i) 
^(^).±(e,p) := Kfi^{e,p) ttnd (ii) fi^^^'^{e,p) := fi^{e,Kp) are spectrally unstable, for all K > Kq. 
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Proof. Note that terms I and III are nonnegative since /ig < 0. So by Lemma 15.21 it suffices for 
instability that the middle integral II be negative and dominate the other two. We begin with 
case (i) of the theorem. Observe that since iJ,^^^'^{e,p) satisfy (|5.6p . the pair {jj,^^^'^ , is 
indeed an homogeneous equilibrium. It suffices to construct a function ij;^ £ iJ^^(D) such that 
(^2^'*! < 0. We choose 



( 1 

'\\}\{r), < r < - 

-1/^1(1 -r), - < r < 1, 



where = 7C(^ ~ ''')) foi^ some normalizing constant 79. Clearly ('0*)/? = so that /// = 0. 

Moreover, 



I 



1 /2 



( y da; > 4coK7r y y (i;)-Vi/(e) dv\il)^\^ rdr 

■oKtt(^ [{vy^vjuie) dv) I r^\i)^\^ dr. 



by choice of the constant 70. So it remains to show that II < —1. 
Using the assumption (|5.8|) and the first scaling (i), we have 

-II = 2 

> AcoKtt 

The integral in f is a finite positive constant thanks to the decay assumption on /x. So we can 
choose K large enough that — // > 1. This settles case (i). 
Similarly, for case (ii) we have 

-11 = 2 {vy\Kp)npie,Kp) dv^lij.l'' dx 

> AcqK^tt j {vy^p^u{e) dv\ilj^\^ rdr = AcoK^tt(^ j {v)-\ju{e) dv) r^lV'*!^ dr, 

which is again greater than one for sufficiently large K. This settles case (ii). □ 

We remark that the constant Kq in Lemma 15.31 is certainly not optimal. For instance, we could 
take V'* to be the ground state of the operator — A^., which is a Bessel function. 



5.2.2 Inhomogeneous equilibria 

For spatially dependent equilibria we will prove a similar result. We first observe as in the homo- 
geneous case that the projection V~ satisfies 

'P'i^) = - [ Hr) dx = 2{ri;)R, (5.10) 
71" Jd 
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for functions ^ = ip{r), where := Jq ip{r) dr. We recall from (|5.2p that 

(>i^^,V')L2 =^(|5,V'|' + ^|V'|') -2^( J rve^tp dv)\^\^ dx-2 j^j {p- {ve^^ dvdx 

= I + II + III, where clearly / > and /// > 0. Of course, both e = {v) and p~ = r{vQ — ip^) 
belong to the kernel of D~. Thus we have 

p-ivetp) = {v)-^V-{ve^) = {v)-^V- (r{ve - + V'V) = {v)-^P~V~ (^) + 
Since Tp^ and are functions depending only on r, we apply (IS.lOp to give 

r^iveij) = 2{vy^p~^R + 2(7;)-i(rV»R. (5.11) 
By definition of p~ , we can write 

rvgfip dv = J {v)^^p^fip dv + rip^ J {v)^^fip dv. 
Using the inequality (a + 6)^ < 2a^ + 26^, we then obtain 

{Ali;,^P)L2< [ + dx-2f n {vy^p-^ip dv)\,p\Ux 



(v) ^(p )Ve dvdxllpRl"^ + 2 sup ( (v) Vp M^^j II "01112(15) 

re[0,l] ^ 

+ 16 sup ( [ {v)-'^\fi-\ dv){r'ip^'ip)l = I + IIA + IIIA + IIB + IIIB. 

■re [0,1] ' 

(5.12) 

We now scale in the variable p to get the following result. 



Theorem 5.4. Assume that ji^ satisfy (|5.6p and i/ia^ 

PlJ-pie-.p) > co/j^(e), V e,p, (5.13) 

for some positive constant cq and some nonnegative function z/(e) such that ^ 0. Define 
H^^^'^{e,p) := iJ,^{e,Kp) and let ip^^^'^ be the solution of the equation 

-A,^W'0 = f ve L(^)'+((f), r(t;e + ^(^)'0)) - ^Ji^^^^~ {{v) ,r{ve - V'^^^'^))] dv, (5.14) 



wzi/i ■(/;(^)''^ = on the boundary 911). T/ien </iere exists a positive number Kq such that the inho- 
mogenous purely magnetic equilibria {/j,^^^'^ ,0, B^^^''^) , with B^^^'^ = ^dr{rip^^^''^), are spectrally 
unstable for all K > Kq. 
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Proof. As before, we will check the instability by showing {A2ip^,'ip^) l2 < for some V*- As in 
Lemma [531 we make the simple choice of ^* given by (j5.9p so that = 0, whence III A = 0, 

with the constant chosen so that the first integral term / in (j5.12p equals 1. Thus we obtain 



+ Coi^||V'(^)'°||i- sup n{vy'\f^~{{v),Kp~)\dv) (5.15) 



rG[0,l] 

+ Co||V^^')'°||i^ sup ( [{v)~^\f,~{{v),Kp~)\dv), 

rG[0,l] ' 

for some constant Cq that depends only on the l? norm of il)^. We shall show that the second integral 
11 A in this estimate dominates all the other terms if K is large. Prom the decay assumption (jl.l2p 
on /i^ and /x^ , we have 



dv<C,j ^^^^^^^^^^^ dv < C„ 

with 7 > 2 and for some constant independent of K. A similar estimate holds for the last 
integral in (jS.lSp . Now by using the assumption ()5.13p and the fact that v{{v)) is even in vq, we 
have 



11A = -2J {v)'^Kp~fi~{{v),Kp-) dv^ dx 

< -2coK^ j ( y r2(^;)-i(t;e-V'(^)'0)2j,((^)) dv^i^^lUx 

= -2coK^ j ( y {v)'\ju{{v)) (i?;)rV*l' dx-2coK^(^ j {v)-^v{{v)) dv) j r^jV^^^)'"! V* 

< -2coA'2(^ / {v)~'^vlu{{v)) dv)\\rip4\2(p) = -ciK^\\rip4\2^^y 



? dx 



where ci > is independent of K. Combining these estimates, we have therefore obtained 

Furthermore, the I?' norm of r-f/^^, is nonzero. We claim that tjj^^^'^ is uniformly bounded indepen- 
dently of K. Indeed, recalling that ip^^^'^ satisfies the simple elliptic equation (j5.14p and using the 
decay assumption (I1.12p on /i^, we have 



|A.^W'0| <C,I dv < C, 



for some constant independent of K. Thus letting u^^\r) := ■0(^)''^(r) + C^r^/S, we observe 
that —ArU^^^ < in D and u^^^ = on the boundary dO. By the standard maximum principle, 
is bounded above and consequently so is same way they are bounded below. 

This proves the claim. Summarizing, we conclude that (^2''/'*5 is dominated for large K by 

11 A and it is therefore strictly negative. □ 
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A Equilibria 



This appendix contains (i) the proof of regularity of our equihbria and (ii) the construction of some 
simple examples of equilibria. We recall that E*^ = —dr^^er and = ^dr{rip^) where {ip^jip'^) 
depends only on r and satisfies the elliptic ODE system (jl.lip . which we rewrite as 



Aip' = h{r, (^0, ^0), -A.V'" = 9{r, V'"), 

({v) + r{ve + V'")) - 1^- {{v) - ^^,r{ve - 



(A.l) 



dv 



V0 



({v) + r{vg + ) - ^l~({v) - T{ve - 



dv. 



For the regularity (i), we will verify that (p^,ip^ G C(B) implies that 'E^,B'^ G C-'^(D). Observe 
that (f^ = 1 and = logr are two independent solutions of the homogeneous ODE Aip^ = 
with wronskian 1/r. Similarly, iIj^ = r and ip^ = 1/r are two solutions of the homogeneous ODE 
A^-i/^^ = with wronskian — 2/r. Thus all the solutions ((/j*^, '0'^) to (|A.ip satisfy the integral 
equations 

nr 

ip^[r) = a+ s{logs — logr)h{-,ip^ ,'il)^){s) ds + ^logr, 
Jo 



/3r + 



2r 



(A.2) 



with arbitrary constants a, /?, 7, 6. Since if^ and are assumed to be continuous at the origin, 
we require 7 = 5 = 0. Clearly h is continuous in D, so that £ C2(D) by ([A^ and E° = 
G C"^(D). As for g, we note that limj,_j.o+ g{r,ip^{r),ip^{r)) = 0, which follows from the fact 
the integrand is odd in vg. So g is also continuous in all of B. Hence, = ^dririjj^) G C^(D), as 
can be seen from (jA.2p . 

As for (ii), the construction of some equilibria, for simplicity we only consider the case ii+{e,p) = 
^"(^iP) = which we take to be an arbitrary function subject to the conditions in (jl.l2p . 

Of course, E'^ = 0, B*^ = is automatically an equilibrium for any fi. However, let us consider 
the inhomogeneous case. Clearly, h{r, 0, 0) = g{r, 0, 0) = 0. No boundary condition is required on 
if^jip^. It is easy to choose so that the functions h{r,-,-) and g{r,-,-) are uniformly bounded, and 
so that for all ^i,?7i,^2!^2 they satisfy 



(A.3) 



for some 9 < 1. Such an assumption is satisfied for instance if /i is uniformly Lipschitz in its 
variables and // replaced by for sufficiently small e. 

Now we denote by T{(p^,tp'^) the right sides of the integral equations in ()A.2p with 7 = 5 = 0. 
It is clear from assumption (jA.Sp that T is well-defined from C([0, 1]) x C([0,1]) into itself. In 
addition, T is a contraction map on a small ball B in this space if a and /3 are sufficiently small. 
So for each small a and (3 there exists a unique solution {^p°,tp°) in B to (TO]) and thus to (|XTj) . 
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